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Abstract. We investigate the influence of a dipole interaction with a classical 
radiation field on a qubit during a continuous change of a control parameter. In 
particular, we explore the non-adiabatic transitions that occur when the qubit is 
swept with linear speed through resonances with the time-dependent interaction. 
Two classical problems come together in this model: the Landau-Zener and 
the Rabi problem. The probability of Landau-Zener transitions now depends 
sensitively on the amplitude, the frequency and the phase of the Rabi interaction. 
The influence of the static phase turns out to be particularly strong, since this 
parameter controls the time-reversal symmetry of the Hamiltonian. In the limits 
of large and small frequencies, analytical results obtained within a rotating-wave 
approximation compare favourably with a numerically exact solution. Some 
physical realizations of the model are discussed, both in microwave optics and 
in magnetic systems. 



PACS numbers: 32.80.Bx, 75.10.Jm, 32.80.Qk, 42.50.Dv 

1. Introduction 

An essential ingredient to a quantum computer is a set of parameters that is 
controllable in the sense that it is possible to manipulate the parameter values at 
any time such that the qubits undergo one-qubit or two-qubit gate operations. For 
quantum computer implementations that rely on nuclear magnetic resonance £Q or 
on spins in quantum dots [2], such a manipulation is possible by switching magnetic 
fields that act on the qubit. This includes the possibility of inverting the sign of an 
acting magnetic field. As the field changes sign, the diabatic energy levels of the 
qubit typically cross. If at the same time a second magnetic field acts in any other 
direction, the adiabatic levels form an avoided crossing instead of an exact crossing. 
Then, depending on the speed at which the control parameters are manipulated, the 
state of the qubit can follow the adiabatic energy levels or undergoes a non-adiabatic, 
so-called Landau-Zener (LZ) transition to the opposite branch |!3 0] 03 EH • 

In the context of quantum computation, it has been proposed to exploit LZ 
transitions for improving the readout of qubits via the so-called Zener flip quantum 
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tunneling This mechanism has recently been implemented for flux qubits HJ. 
A method for non-adiabatic electron manipulation in quantum dots also relies on 
LZ transitions (TJ3- Moreover, the observation of LZ transitions is a clear sign of 
coherence like, e.g., optical coherence in a classical optical ring resonator JT] or 
macroscopic quantum coherence in superconducting loops ^l^H- LZ transitions have 
also been used to determine tiny interactions between levels in molecular clusters |14| . 
While in these cases, LZ transitions are beneficial, the opposite is true in the case 
of adiabatic quantum computing ^JE|- There, the computation is performed by a 
quantum system that follows adiabatically the instantaneous ground state of a slowly 
varying Hamiltonian and, consequently, the emergence of any non-adiabatic transition 
constitutes an error source. 

The physical origin of a coupling between two levels of a quantum system is not 
necessarily simply an overlap between the respective wave functions. In particular for 
spins and atoms, such a coupling typically stems from the dipole interaction of the 
system with a radiation field. In a seminal work J7| (see also Rabi predicted 

within an exact quantum mechanical treatment that a classical, monochromatic and 
circularly polarized radiation field induces spin rotations with a frequency which at 
resonance is proportional to the field amplitude. With resonant linearly polarized light, 
the same characteristic harmonic Rabi oscillations of atomic inversions are observed. 
As a consequence of the linearly polarized driving, the optical realization of the Rabi 
problem is not exactly solvable. For resonant excitations, however, it is possible to 
apply a rotating-wave approximation (RWA) which formally restores the situation 
with circular polarization I2()| . This necessarily neglects effects beyond RWA like 
the Bloch-Siegert shift of the resonance frequency. In optical realizations, however, 
such these effects are very tiny ^] . 

The question now arises whether a level interaction mediated by a classical 
monochromatic radiation field can induce Landau-Zener transitions in a two-level 
system as its energies cross. In this work, we demonstrate that this is indeed the 
case. Thereby, we investigate LZ transitions that are induced by the coupling of a 
spin to a linearly polarized light field, henceforth referred to as Rabi coupling. In the 
traditional LZ problem, non-adiabatic transitions occur when the adiabatic energy 
levels are close to each other. By contrast, we will find that with a Rabi coupling to a 
high-frequency field, the transitions take place at times at which the radiation field is 
at resonance with the diabatic energy levels. This allows for sufficiently weak coupling 
a perturbative treatment within a rotating-wave approximation. For suitably chosen 
parameters, the driving reduces the probability for LZ transitions which relates this 
problem to the so-called coherent destruction of tunneling (23 EH ■ A different kind of 
time-dependence would be provided by coupling the two-level system to a noise source. 
In this paper, we will work in the coherent limit instead, which at least in possible 
realizations of quantum computation should be a good approximation on a relevant 
time scale. In relation to this, it is interesting to note that LZ tunneling is fairly robust 
against classical noise HUES EH EH and quantum dissipation [2lfl EM EM EH1 1521 1551153] . 

In this context, we like to emphasize that herein considered Landau-Zener 
tunneling due to coupling to a light field is different from the one considered in 
Refs. |24l 1851 ET5] where it is the diabatic energies of the two-level system that are 
subject to a time-periodic modulation (the so-called dynamic Stark effect). Still, the 
model considered here applies quite generally to qubits whose level interaction varies 
harmonically in time due to interaction with an external field that causes a negligible 
dynamic Stark effect. In the discussion at the end of the paper some possible physical 
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realizations of the model are suggested. 



2. The Landau-Zener model with harmonic interaction modulation 

We consider a quantum system ("atom") with two relevant energy levels |1) and |2) 
whose time-dependent energies ±Vt/2 cross at t = 0. Both levels are coupled to 
a classical dipole field with frequency f2 and phase <f>. The effective amplitude g is 
given by field strength times the dipole moment of the two-level system. Thus, the 
Hamiltonian reads 

H{t) = ^a z + f{t)a x , /(t)=0cos(fif + 0, (1) 

where cr z \l) = |1) and cr z \2) = — 12). Moreover, we assume that at time t = — oo, the 
system is in its instantaneous ground state, i.e. \ip(— oo)) = |1). If the energies in the 
first term would be time- independent, and if the field would be in resonance with the 
atomic energy difference, then this Hamiltonian would describe an atom undergoing 
Rabi oscillations. 

Since for most times, the Hamiltonian is dominated by its first term, a 
proper interaction- picture representation is provided by the transformation Uo(t) — 
ex.p(-iVt 2 <r z /2h), that is \tp(t)) = U {t)\ip(t)} and \4>(t)) = [ci(t)|l> +c 2 (t)\2)], where 
the interaction-picture probability amplitudes obey 



Cl 



f{t) e iVt2 / 2h 



c 2 J h\ f(t) e~ iVt ' 2h 




(2) 



For fi = cf> = 0, the Hamiltonian defines the standard Landau-Zener problem 
for which the exact solution of the equation of motion (J2J) can be expressed in terms 
of parabolic cylinder functions [3]. Then, the time-evolution from t — — oo to i = oo 
is given by the S-matrix 



where q = exp(—2irg 2 /hV) and the Stokes phase \ = 7r /4 + argT(l — id) +6(lnS — 1), 
with 5 = g 2 /(hV) and T(..) the Gamma function. The famous Landau-Zener transition 
probability follows readily: the probability P that the atom ends up in the initially 
unoccupied level |2) is 

P= |c 2 (t = oo)| 2 = 1 - e" 2 ^ 2 ^. (4) 

Note that this result is exact for all values of g and V. 

Already a good approximation to the time-dependent solution is provided by the 
fact that for f(t) = g and sufficiently large times, the phase factors in the matrix in 
Eq. J2Jl are rapidly oscillating with a quadratic time dependence. As a consequence, 
ci and £2 remain essentially constant. By contrast, close to t = these phase factors 
assume a stationary value and the two-level system undergoes a transition. This means 
that S-matrix @ in fact describes a transition taking place at t = 0. Thus, in the 
interaction picture the dynamics is approximately given by \ip(t)} = \ip(— oo)) for t < 
and \$(t)) = S g |-0(oo)) for t > 0. 

For 7^ 0, the time t = no longer marks the time at which the phase is 
stationary and the behaviour changes significantly, as we will see below. In order to 
anticipate the richness of the resulting dynamics, we have numerically integrated the 
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Figure 1. LZ transition probability |c2(t)| 2 as a function of time in units of 
(h/V) 1 ' 2 for <f> = and various values of the interaction strength g in units 
(HV) 1 / 2 and the modulation frequency Q in units of (V/h) 1 / 2 . The red lines 
mark the standard LZ transition probability [Eq. J1J, dashed], the RWA result 
for a double transition [Eq. solid], and the transition probability at the 

intermediate stage [Eq. I1UI . dotted]. The green curves in the four upper-left 
panels correspond to the adiabatic-following result Eq. 



equations of motion (J2J for (f> = and various coupling strengths g and frequencies O. 
Fig- depicts the time dependence of the probability |c 2 (i)| 2 to find the atom in the 
initially unpopulated level |2). For small interaction g (see lower five plots in Fig.^J, 
the final transition probability for long times does not depend strongly on frequency, 
although the curves differ strongly around t — 0. The most interesting feature of 
Fig. ^ is that for high frequencies, ^> \/V/h, the dynamics consists of two (almost) 
independent transitions. 



3. Adiabatic vs. non-adiabatic transitions 

In the standard Landau-Zener problem, one distinguishes two limiting cases: If the 
level crossing occurs very rapidly, the potential switches practically instantaneously 
such that no significant dynamics can take place. The system will then remain in level 
|1) so that finally P = 0. In the opposite limit the instantaneous energy levels change 
very slowly. The system then follows adiabatically the lower-energy level \E-(t)) and 
ends up with P = 1. 

3.1. Adiabatic following 

Adiabatic following means that transitions between the instantaneous eigenstates 
can be neglected. The criterion for adiabatic following is that at each instance of 
time, the coupling between the adiabatic energies is "sufficiently small" , much smaller 
than the energy splitting. Stated in mathematical terms, this requirement becomes 
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Figure 2. Adiabatic energies E±(t) for a Rabi coupling with <f> = (solid lines) 
and (f> = 7r/2 (dashed). The dotted lines marked the limiting values for g = and 
the standard LZ model with constant coupling, respectively. The other parameters 
are g = 4vW, and Q = 0.5-fV/h. 



|( E-(t)\^\E + {t) )| < - S_(t)|/fi. This gives for the standard Landau-Zener 

problem a splitting 2g and the adiabaticity condition HV <C <7 2 . 

For the time-dependent two-state Hamiltonian (JJJ, the condition for adiabatic 
following becomes more involved, because the minimal splitting of the adiabatic 
energies E±(t) = ±y/ (Vt/2) 2 + f 2 (t) depends not only on the coupling strength g 
but also on the frequency SI and the phase 4>. The minimal splitting will never be 
larger than 2g, whatever the frequency and phase. The sensitive dependence on the 
phase becomes obvious from Fig. El In particular for <f> = ir/2, the adiabatic spectrum 
no longer exhibits an avoided crossing but rather an exact crossing, since at time 
t = both terms in the Hamiltonian Q vanish simultaneously. Consequently, the 
adiabaticity condition is violated irrespective of the values of the g and f2. This 
qualitative difference already provides a hint that the phase (j) has a strong influence 
on the population dynamics. 

By contrast, for a phase </> ^= 7r/2(mod7r), the energies E± never form an exact 
crossing. Thus, it is always possible to choose V and fl so small that the adiabaticity 
condition is fulfilled. Then, by making the time-dependent transformation to the 
instantaneous-energy representation, neglecting there the off-diagonal elements in the 
equations of motion, solving the dynamics and then transforming back to the diabatic 
representation, it follows that the probability | C2 ) | 2 to find the qubit in diabatic state 
1 2) goes from zero to one as 

\c2,a.diabaUc{t)\ 2 = ^ ( H j = = j ■ (5) 

V ^{Vtf +Ap{t)) 

The dependence for intermediate times on the interaction modulation f(t) is clearly 
seen. The four upper-left panels in Fig. [I] with £1 = or 0.5 and g — 1 or 4 are very 
well described by the adiabatic- following result JSJ. 
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3.2. Non-adiabatic regime 

For phase = 0, the Rabi coupling f(t) is zero at time t = tt/20, and the energy 
splitting becomes nV/2$l. This means that for a large driving frequency > yJV/h, 
the adiabaticity condition is violated. The data shown in the right columns of Fig. ^ 
indicate that in this regime the dynamics consists of two transitions at times ^Ml/V. 
If the time 2MI/V between the individual transitions is sufficiently large, as specified 
below, the two transitions are essentially independent of each other. Then, it is possible 
to derive within a rotating-wave approximation an analytical expression for the final 
transition probability. The derivation is closely related to the transfer matrix method 
employed in Refs. (571 l2l) . 

With the new variables di{t) = c ± cxpiiMl 2 / AV) and d 2 it) — c 2 exp(— ihVl 2 /W), 
one obtains from Eq. @ the equations of motion 

9 



dl= - ' l 2h 

a ■ 9 

di = — l — 

2h 



e iv(t+hn/vy /2h , e iv(t-nn/vy/2h 



-iv(t+nn/v) 2 /2h _^ e -iv(t-nn/v) 2 /2fi 



h, (6) 
di. (7) 



Like in Eq. J2J, the phases on the right-hand side obey a quadratic time dependence. 
Thus, with the arguments provided after Eq. Q, we can conclude that each phase 
factor is relevant only at times at which the phase is stationary, i.e., the first term 
contributes only at time i_ — —Ml/V while the second term becomes relevant at time 
t_l_ = Ml/V. Thus, we keep at both times t_ and t+ only the respective resonant term 
while the "counter-rotating" term is neglected. (To be sure, two separate rotating- 
wave approximations are needed corresponding to t + and t_.) Then, at times close 
to t^, the equation of motion is of the same form as Eq. J3J) and the dynamics is 
determined by the S-matrix J2J) with g replaced by g/2, i.e. S T = S g / 2 - Consequently 
within the transfer matrix approximation, the time evolution becomes 

( |^(-oo)) for t < -Ml/V, 

\{j 1 (t)) = \ S 3/2 |^(-oo)) for -Ml/V < t < Ml/V, (8) 

{ Sg/aM-oo)) for t > Ml/V. 

With this expression, the probability to find the system at time t = oo in state \2) is 
readily evaluated as 

P = |<2|S" /2 |1)| 2 = 4 e -^ 2 / 2HV (l - e~^ 2 / 2hV ). (9) 

During the intermediate times — Ml/V < t < Ml/V, the occupation probability of 
level 1 2) becomes 

Pint - l(2|S 9/2 |l)| 2 = 1 - exp(-7r. g 2 /2r^). (10) 

Note that in these expressions, the exponent differs from the exponent in Eq. 
by a factor 1/4. Moreover, P in Eq. @ no longer depends monotonously on the 
coupling strength as in the standard LZ problem, but rather assumes a maximum 
for exp(— Trg 2 /2hV) = \. Interestingly enough, the transfer matrix results are 
independent of the Stokes phase \ and the modulation frequency tt. The independence 
of the frequency is confirmed by Fig. Q where the transfer-matrix results (red lines) 
nicely agree with the exact results (blue lines) for SI — 2, 10, and 40. Clearly, at long 
times, the probability |c 2 | 2 is a function of only the coupling strength g. 

Below we will compare more systematically the transition probabilities obtained 
from the transfer matrix method with a numerically exact solution. But first we have 
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Figure 3. Final transition probability P as a function of coupling g, for several 
values of the interaction modulation frequency Q. The dashed line marks the 
standard LZ transition probability J1J valid for Q = 0. The transfer matrix 
result coincides with the numerical result for Q = 40. (The numerical time 
integration was performed from times — 500i/ h/V to +500\/ fo/V.) 

to specify the conditions under which the time between two LZ transitions will be 
long enough for the transfer-matrix analysis to hold. It should be remembered that LZ 
transitions neither occur instantaneously nor take infinitely long [331 [3011121- The time 
2MI/V between the two consecutive transitions should be larger than the duration of 
a single transition. It has been estimated that the standard LZ transition has a typical 
duration tlz — yJh/V when non-adiabatic transitions are probable (2g 2 /(hV) <C 1), 
while tlz — 2g/V in the adiabatic limit 2g 2 /(hV) 1 |38| ■ Correspondingly, reliable 
results of the transfer matrix approach are to be expected if 

n £ y/V/(4h) for g 2 /(2hV) < 1, (11) 

0£s/ft f or g 2 /(2hV)^> 1. (12) 

These estimates are confirmed by the numerical solution of Eq. J2J as plotted in Fig. [31 
which depicts the probability of finding the system in state |2) at large times. The 
figure makes clear that for small coupling strengths g <C \/KV the condition is 
sufficient but not necessary, because Eq. © is seen to be accurate irrespective of the 
frequency. This is in accordance with the lower five panels in Fig.^ an d with fact that 
the standard LZ result and the transfer-matrix expression @ agree that P equals 
2irg 2 /V to first order in g 2 /V. On the other hand, we find significant deviations from 
the expression © once g becomes larger than y/hV and of the order hfl. This is where 
the two LZ transitions start to "feel" each other. For sufficiently large coupling, the 
transition probability even increases again and assumes further maxima with P = 1. It 
is nice that the argument can be turned around and that LZ times can be estimated 
with the help of the frequencies at which the exact and the transfer-matrix results 
start to deviate. In doing so, we indeed find (here and in Sec.^J that tlz ~ 2g/V for 
2g 2 /(hV) 3> 1, in agreement with [3E]- The present model provides an independent 
and simple method to determine Landau-Zener times. 
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Figure 4. Final transition probability P as a function of coupling strength g and 
frequency Q. Blue areas correspond to P = 1, white areas to P = 0. 

4. Exploring the crossover region 

Our analysis has identified two different parameter regimes in which the analytical 
solution is confirmed by the our numerical results. First, there is the regime of slow 
driving in which Ml denotes the smallest energy scale of the problem. Then, the 
time dependence of the coupling is not essential and the transition probabilities are 
the same as in the standard LZ problem. In the second regime Ml is the largest 
energy scale and the transfer matrix results hold. In particular, we find P = 1 for 
irg 2 = 2KV\u2. In this section, we complement our analytical findings by numerical 
results for the intermediate parameter regime. 

Figure 21 shows the final transition probability, that is the occupation of state |2) 
in the limit t — > oo. (The time interval for numerical integration was chosen the same 
as for Fig. |3J) The vertical stripe with P = 1 for fl < 0.5 corresponds to the adiabatic 
regime. The horizontal blue stripe marks the maximum found within the transfer 
matrix approach. The figure also confirms (i) that the location of the maximum has 
no significant frequency dependence and (ii) that P decays for a larger coupling g 
almost to zero, yielding the white region with P = above the horizontal blue band, 
in agreement with the transfer-matrix prediction JjJJ. Increasing g further, we find 
that at g « Ml, the transition probability again assumes values close to unity. This 
regime, including the sequence of maxima and minima that can be observed for even 
larger coupling, is beyond the range of validity of the transfer matrix method. The 
fact that the transfer matrix approach starts to break down along the diagonal g = Ml 
in Figure 01 neatly agrees with the estimate l|12l) . With the reasoning given in Sec. 13.21 
we can infer from Fig. 21 that the estimate for the LZ time tlz ~ 2g/V holds at 
least in the broad parameter regime 1 < g 2 /hV < 100. 
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5. Phase dependence 

When discussing the adiabatic energies of the Hamiltonian 0), we have already 
anticipated that the phase cj) might have some relevance that we explore in the 
following. For that purpose, we adapt the analytical approach of Sec . 13 . 21 accordingly. 



5.1. Transfer matrix approach 

Inserting again the definitions di — C\ exp(iM7 2 /4I/) and d 2 = c 2 exp (ifiil 2 /4V) into 
the equation of motion J3J), we find 
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dl= - Sn 
d 2 =-i^- 

2h 



jV(t+Hn/vy /2h+i<t> i e iV(t-hn/vy /2H-i<l> 



e -iV(t+hn/vy /2h-\(j> , e -iV(t-hn/vy /2h+i4> 



(13) 
dx. (14) 



These equations differ from Eq. merely by the phase <fi in the exponents. The goal 
is now to transform Eqs. (| 1 31) and 114|) such that they assume at times t T = ^phfl/V 
the same form as Eqs. and Q). After such a transformation the transfer matrix 
method could be used again. At time t_ , when only the first term in the equations of 
motion is relevant, an appropriate transformation reads 



T = C */» J " < 15 > 
The corresponding transfer matrix S_ follows from a transformation of S g / 2 with T and 
reads S- = T 1 S 9 / 2 T. With the same reasoning, we find that at time t+ the required 
transformation is 1 1 and the transfer matrix is S+ = TS 9 / 2 1 . Consequently, the 
complete time evolution becomes |^(oo)) = S+S_|V'(— oo)). This leads to the final 
transition probability 

P = 4 e"^ 2 /2V / 1 _ e -2^/2v\ CQS 2 ^ (16) 



i<t>/2 Q 



Clearly, the phase shift modifies the transition probability by a factor cos 2 (p. At 
intermediate times t « 0, the occupation is determined by S_. Interestingly enough, 
the absolute values of its matrix elements do not depend on <f> and, thus, we find at 
t = the same phase-independent transition probability P m t as in JTDJ|. Evidently, the 
phase dependence in (|16|l is caused by quantum interference between the two transition 
paths from |1) to |2), which also explains the absence of any phase dependence after 
the first transition. 

Figure [3 shows a comparison with the numerically computed time evolution. The 
long-time limits compare favourably with our prediction for the final state. As for 
the special case = 0, the final occupation is independent of the frequency, provided 
that Q, exceeds g/h and \JV/h. The strong phase dependence in (|16|l may come 
as a surprise, since it is tempting to argue that for high frequencies, phase relations 
should be immaterial due to the many oscillations occurring during each LZ transition. 
However, figure |S] clearly demonstrates that such reasoning is incorrect. 



5.2. Time-reversal anti- symmetry 

The transition probability obviously vanishes for <p — ir/2. This behaviour can 
already been obtained from symmetry arguments. For this phase, the Hamiltonian 
Q is anti-symmetric under time reversion t — > —t, i.e., H(t) = —H(—t). Then 



Landau-Zener transitions in qubits controlled by electromagnetic fields 10 
1 

0.8 
0.6 
I 0.4 
0.2 


-40 -20 20 40 

Figure 5. Transition probability as a function of time for g/y/hV = 
[2 In (2)/7r] 1//2 0.664 and driving frequency £7 = IQ^/V/h. The transfer matrix 
method predicts the intermediate population |c2 (0) | 2 = 0.5 and the final transition 
probabilities 1, 0.75, 0.25, and 0, respectively. 

the time evolution operators U(t, 0) and U{— t, 0) obey the same equation of motion. 
Moreover, they obviously are identical and equal to 1 at t = 0. The equalities 
f7(oo,0) = U{— oo,0) = W{Q,—oo) follow immediately, the last one from unitarity. 
Consequently, we find 

U{oo, -oo) = J7(oo, 0)17(0, -oo) = U+(0, -oo)Z7(0, -oo) = 1, (17) 

which implies that at long times, the system will evolve back to its initial state. 

This ideal back-evolution relates our problem to the Loschmidt echo which has 
been employed for testing the sensitivity of a "chaotic" quantum system on weak 
perturbations |41|. In the present case, the small parameter is 5<f) — (j> — ir/2 
which corresponds to the perturbation Hamiltonian £ = — 2gcos(f2i) sm(6<p)tr x . 
Note however, that the present system does not exhibit any sensitive exponential 
dependence on the perturbation. 

6. Discussion and summary 

Our study of diabatic level crossing in a system subject to a time-dependent dipolc 
force revealed two intriguing features. First, the probability for non-adiabatic 
transitions is not simply a monotonous function of the coupling strength but exhibits 
several maxima and minima. In particular, it vanishes for zero coupling and equals 
unity if the relation gj\fW = [21n(2)/7r] 1 / 2 w 0.664 is fulfilled. This is in contrast 
to the standard Landau-Zener problem where the extreme cases require a vanishing 
or an infinite interaction strength. Second, we found that the phase ef> of the dipole 
field has a significant influence on the transition probability which is proportional to 
cos 2 <j). The combination of both effects enables one to steer the system towards the 
one or the other final state. In turn, it is also possible to use the setup as a diagnostic 
tool for an unknown phase of a radiation field. The fact that the results of the transfer 
matrix approach are only valid if the duration tlz of a single Landau-Zener transition 
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is sufficiently small, provides a further application. Measuring the frequency at which 
the approximation breaks down, allows one to determine tlz . 

A straightforward physical realization of our set-up is naturally provided by spin-i 
systems in time-dependent magnetic fields. Moreover, one could think of experiments 
with effective low-spin systems such as molecular complexes Indeed, the recent 

resonant-photon absorption experiments on the effectively spin-^ molecular complex 
V15 could also be performed with a sweep of the magnetic field that is called B$. A 
further possible realization is given by Rydberg atoms in the vicinity of a crossing of 
the highest Stark level in the nth manifold of the atom and the lowest Stark level of the 
n+ 1st manifold 44, 45 . In this case, the energies are swept by the dc Stark effect and 
the interaction is driven by a harmonically varying microwave field. Then, the strong 
suppression of non-adiabatic transitions can be tested experimentally. Particularly 
promising in this respect would be variants of microwave ionization experiments of 
Rydberg atoms that are based on a mechanism of multiple Landau-Zener transitions 
to higher and higher Stark manifolds |45| . 

To summarize, we have studied Landau-Zener transitions in a two-level atom 
subject to a harmonically time-dependent, Rabi-like interaction. As a main difference 
to the standard Landau-Zener problem, we find that by tuning the coupling strength 
within a relatively small range, it is possible to continuously change the transition 
probability from zero to unity. This behaviour can be explained within a transfer 
matrix approach, which provides reliable results provided that the driving frequency 
is sufficiently large. Moreover, this analytical approach allows one to determine the 
influence of the phase relation between the diabatic energy crossing and the dipolc 
field. It revealed that the transition probability is proportional to cos 2 and therefore 
will vanish for (f> = n/2. The latter result was also shown by analyzing the underlying 
time-reversal symmetry. The sensitive phase dependence can be exploited both for 
steering the system towards a particular state and for measuring an a priori unknown 
phase relation. 
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